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Abstract This paper reports on the microscopic calcu- 
lation of ground and excited states Gamow- Teller (GT) 
strength distributions, both in the electron capture and 
electron decay direction, for 54,55,56^^^ rpj^^ associated 
electron and positron capture rates for these isotopes 
of iron are also calculated in stellar matter. These cal- 
culations were recently introduced and this paper is a 
follow-up which discusses in detail the GT strength dis- 
tributions and stellar capture rates of key iron isotopes. 
The calculations arc performed within the framework 
of the proton-neutron quasiparticle random phase ap- 
proximation (pn-QRPA) theory. The pn-QRPA the- 
ory allows a microscopic state-by-state calculation of 
GT strength functions and stellar capture rates which 
greatly increases the reliability of the results. For the 
first time experimental deformation of nuclei are taken 
into account. In the core of massive stars isotopes of 
iron, 54,55,56pg^ g^j.g considered to be key players in de- 
creasing the electron-to-baryon ratio (Ye) mainly via 
electron capture on these nuclide. The structure of 
the presupernova star is altered both by the changes 
in Ye and the entropy of the core material. Results are 
encouraging and are compared against measurements 
(where possible) and other calculations. The calculated 
electron capture rates are in overall good agreement 
with the shell model results. During the presupernova 
evolution of massive stars, from oxygen shell burning 
stages till around end of convective core silicon burning, 
the calculated electron capture rates on ^^Fe are around 
three times bigger than the corresponding shell model 
rates. The calculated positron capture rates, however, 
are suppressed by two to five orders of magnitude. 
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1 Introduction 

Supernovae of Type la and Type II are considered to 
be the two major contributors to the production of el- 
ements in the universe (see for example, (jll)). Stars in 
mass range M < 8Mq [Mq denotes the solar mass) end 
their life as white dwarfs. The ultimate destiny of more 
massive stars is even more interesting. Once the mass 
of the iron core exceeds the Chandrasekhar limit the 
Pauli principle applied to electrons cannot prevent fur- 
ther gravitational collapse to an even more exotic and 
denser form of matter than in a white dwarf star. A 
violent collapse of the core is initiated which ultimately 
leads to a spectacular supernova explosion (where the 
luminosity of the star becomes comparable to that of 
an entire galaxy containing around 10"'^"'^ stars!). The 
supernova leaves behind a compressed ball of hot neu- 
trons - a neutron star. As the hot neutron star cools, 
any further collapse is prevented by the Pauli principle 
applied to the neutrons, unless the mass is so great that 
the star can become a black hole. The structure of the 
progenitor star has a vital role to play in the mechanism 
of the explosion. Core-collapse simulators, to-date, find 
it challenging to convert the collapse into a successful 
explosion. A lot many physical inputs arc required at 
the beginning of each stage of the entire simulation pro- 
cess including but not limited to collapse of the core, 
formation, stalling and revival of the shock wave and 
shock propagation. It is highly desirable to calculate 
the presupernova stellar structure with the most reli- 
able physical data and inputs. 

A smaller precollapse iron core mass and a lower en- 
tropy should favor an explosion. A smaller iron core 



size implies less energy loss by the shock in photodis- 
integrating the iron nuclei in the overlying onion-like 
structure whereas a lower entropy environment can as- 
sist to achieve higher densities for the ensuing collapse 
generating a stronger bounce and in turn forming a 
more energetic shock wave (l2|). A smaller entropy also 
helps in reducing the abundance of free protons thereby 
lowering the electron capture rates on these free protons 
and resulting in a much higher value of Ye at the time 
of bounce (at this instant of time the iron core collapses 
to supranuclear density and rebounds creating a shock 
wave). As a result, the shock energy 
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is larger for larger final lepton fraction Ye/ (prior to 
collapse) and larger difference of initial and final lepton 
fraction. 

Electron/positron captures and /3 -decays are amongst| 
the most important nuclear physics inputs that deter- 
mine both the Yef and the entropy at the presupernova 
stage. These nuclear weak-interaction mediated reac- 
tions directly affect the lepton-to-baryon ratio. Further 
the neutrinos and antineutrinos produced as a result 
of these nuclear reactions are transparent to the stellar 
matter at presupernova densities and therefore assist in 
cooling the core to a lower entropy state. These weak- 
interaction rates are required not only in the accurate 
determination of the structure of the stellar core but 
also bear significance in (explosive) nucleosynthesis and 
element abundance calculations. As such it is impera- 
tive to follow the evolution of the presupernova collapse 
stage with a sufficiently detailed and reliable (micro- 
scopically calculated taking into account the nuclear 
structure details of the individual nuclei) nuclear reac- 
tion network that include these weak-interaction medi- 
ated rates. Weak interactions in presupernova stars are 
known to be dominated by allowed Fermi and Gamow- 
Teller (GT) transitions (jj). The exact r-process yields 
from Type II SNc are not known and arc related to 
the complete understanding of the mechanism of these 
supernova explosions. Electron capture process is one 
of the essential ingredients involved in the complex dy- 
namics of core-collapse supernovae and a reliable esti- 
mate of these rates can certainly contribute in a better 
understanding of the explosion mechanism. 

SNe la are thought to be the explosions of white 
dwarfs that accrete matter from binary companions and 
are intensively investigated for their accurate calibra- 
tion as the cosmological standard candles (one needs to 
predict the SN la peak luminosity at least with a 10% 
precision) . These explosions are also important for stel- 
lar nucleosynthesis and are strongly contested to be the 
main producers of Fe peak elements in the Galaxy (jsl). 



The abundance of the Fe group, in particular of neutron 
rich species, is highly sensitive to the electron captures 
taking place in the central layers of SNe la explosions 
(jsl). These captures drive the matter to larger neutron 
excesses and thus strongly influence the composition 
of ejected matter and the dynamics of the explosion. 
It is highly probable that electron capture occurs in 
the burning front driving the matter to large neutron 
excess. In particular, GT properties of nuclei in the 
region of medium masses around A=56 are of special 
importance because they are the main constituents of 
the stellar core in presupernova conditions. 

The GT transition is one of the most important nu- 
clear weak processes of the spin-isospin (err) type. This 
charge exchange transition is not only subject of inter- 
est in nuclear physics but also of great importance in 
astrophysics. The GT transitions are involved in many 
weak processes occurring in nuclei, e.g. nucleosynthe- 
sis and stellar core collapse of massive stars preced- 
ing the supernova explosion. In the early stages of the 
collapse, the electron capture and /3-decay are impor- 
tant processes. These reactions are dominated by GT 
(and also by Fermi) transitions. Electron (positron) 
capture rates are very sensitive to the distribution of 
the GT+ (GT_) strength. In the GT+ strength, a pro- 
ton is changed into a neutron (the plus sign is for the 
isospin raising operator (i+), present in the GT ma- 
trix elements, which converts a proton into a neutron). 
On the other hand the GT_ strength is responsible for 
transforming a neutron into a proton (the minus sign 
is for the isospin lowering operator (t_) which converts 
a neutron into a proton). The total GT+ strength is 
proportional to the electron capture strength (jj). 

GT+ strength distributions on nuclei in the mass 
range A = 50 - 65 have been studied experimentally 
mainly via (iijp)charge exchange reactions at forward 
angles (e.g. (|i;li;@)). Similarly GT_ strength distri- 
butions were studied keenly using the (p,n) reactions 
(e.g. P; lZI; Isl)). It has been shown that for (p,n) and 
(n,p) reactions the 0° cross sections for such transitions 
are proportional to the squares of the matrix elements 
for the GT /3 decay between the same states (e.g. ([9|)). 
Results of these measurements show that, in contrast to 
the independent particle model, the total GT strength 
is quenched and fragmented over many flnal states in 
the daughter nucleus. Both these effects arc caused by 
the residual interaction among the valence nucleons and 
an accurate description of these correlations is essential 
for a reliable evaluation of the stellar weak interaction 
rates due to the strong phase space energy dependence, 
particularly of the stellar electron capture rates. 

Fuller, Fowler and Newman (FFN) (jlOl ) performed 
the first-ever extensive calculation of stellar weak rates 



including the capture rates, neutrino energy loss rates 
and decay rates for a wide density and temperature do- 
main. They performed this detailed calculations for 226 
nuclei in the mass range 21 < ^ < 60. They stressed 
on the importance of the Gamow- Teller (GT) giant res- 
onance strength in the capture of the electron and esti- 
mated the GT centroids using zeroth-order {Ohoj ) shell 
model. Both the decay and capture rates are expo- 
nentially sensitive to the location of the GT centroid 
(111). The location of the GT resonance affects the stel- 
lar rates exponentially, while the total strength affects 
them linearly (jllr i. Few years later Aufderheide et al. 
(J12f ) extended the FFN work for heavier nuclei with A > 
60 and took into consideration the quenching of the GT 
strength neglected by FFN. They tabulated the 90 top 
electron capture nuclei averaged throughout the stellar 
trajectory for 0.40 < Ye < 0.5 (see Table. 25 therein). 
Later the experimental results of Refs. (|j; |5|; M, lZ|; IsJ) re- 
vealed the misplacement of the GT centroid adopted in 
the parameterizations of Ref. (jlOl ). Since then theoreti- 
cal efforts were concentrated on the microscopic calcula- 
tions of GT strength distributions and associated weak- 
interaction mediated rates specially for iron-regime nu- 
clide. Large-scale shell model (LSSM)(e.g. (|l3[)) and 
the proton-neutron quasiparticlc random phase approx- 
imation theory (pn-QRPA) (e.g. (jlj)) were used ex- 
tensively and with relative success for the microscopic 
calculation of stellar weak rates. Monte Carlo shell- 
model is an alternative to the diagonalization method 
and allows calculation of nuclear properties as thermal 
averages (e.g. (Il5h). However it does not allow for de- 
tailed nuclear spectroscopy and has some restrictions in 
its applications for odd-odd and odd-A nuclei. 

The pn-QRPA theory is an efficient way to gener- 
ate GT strength distributions. These strength distri- 
butions constitute a primary and nontrivial contribu- 
tion to the capture rates among iron-regime nuclide. 
The usual RPA was formulated for excitations in the 
same nucleus. Halblcib and Sorenson (|l6l ) generalized 
this model to describe charge-changing transitions of 
the type (Z, iV) ^ (Z ± 1, A^ =F 1) and pn-QRPA first 
came into existence more than 40 years ago. The model 
was extended to deformed nuclei (using Nilsson-model 
wave functions) by Krumlindc and Moller (jlTJ) . Exten- 
sion of the model to treat odd-odd nuclei and transi- 
tions from nuclear excited states was done by Muto and 
collaborators (jlSj). 

Nabi and Klapdor-Kleingrothaus used the pn-QRPA 
theory to calculate the stellar weak interaction rates 
over a wide range of temperature and density scale for 
sd- (jlQl ) and fp/fpg-shell nuclei (|l4| ) in stellar matter 
(see also Ref. (l20l)). These included the weak interac- 
tion rates for nuclei with A = 40 to 44 (not yet calcu- 
lated by shell model). Since then these calculations 



were further refined with use of more efficient algo- 
rithms, computing power, incorporation of latest data 
from mass compilations and experimental values, and 
fine-tu ning of model parameters (2lJ; [2^ |23|; |2J; |25|; 
26l : I27I : l28t I29I : l30f ). There is a considerable amount of 
uncertainty involved in all types of calculations of stel- 
lar weak interactions. The uncertainty associated with 
the microscopic calculation of the pn-QRPA model was 
discussed in detail in Ref. (|27l ) . The reliability of the 
pn-QRPA calculations was discussed in length by Nabi 
and Klapdor-Kleingrothaus (|l4l ). There the authors 
compared the measured data (half lives and B(GT±) 
strength) of thousands of nuclide with the pn-QRPA 
calculations and got fairly good comparison. 

Three key isotopes of iron, 54,55,56pg^ were selected 
for the calculation of GT± strength distributions and 
associated stellar electron and positron capture rates 
in this phase of the project. Whereas sufficient ex- 
perimental data are available for the even-even 54,56pg 
isotopes to test the model, ^^Fe has low-lying excited 
states which have a finite probability of occupation in 
stellar conditions and a microscopic calculation of GT-t 
strengths from these excited states is desirable. Aufder- 
heide and collaborators (|l2f ) ranked 54,55,56pg amongst 
the most influential nuclei with respect to their impor- 
tance for the electron capture process for the early pre- 
supernova collapse. Later Hcger ct al. (|3lr ) studied 
the presupernova evolution of massive stars (of masses 
15M0, 25M0, and 40Afo) and rated 54,55,56pg amongst 
top seven nuclei considered to be most important for 
decreasing Y^ in 15Mq and 4OM0 stars. In 25M0 stars 
these isotopes of iron were ranked as the top three key 
nuclei that play the biggest role in decreasing Y^. These 
isotopes of iron are mainly responsible for decreasing 
the electron-to-baryon ratio during the oxygen and sili- 
con burning phases. Besides, ^^Fe was also found to be 
in the top five list of nuclei that increase Ye via positron 
capture and electron decay during the silicon burning 
phases. This paper presents the detailed analysis of the 
improved microscopic calculation of GT-t strength dis- 
tributions. Details of calculation of stellar electron and 
positron capture rates for these three isotopes of iron 
using the pn-QRPA model are also furnished in this 
manuscript. Comparisons against previous calculations 
are also presented. These improved calculations were 
introduced recently in Ref. (1281 ) where it was reported 
that the betterment resulted mainly from the incorpo- 
ration of measured deformation values for these nuclei. 
The idea is to present an alternate microscopic and ac- 
curate estimate of weak interaction mediated rates to 
the collapse simulators which may be used as a reliable 
source of nuclear physics input to the simulation codes. 

The next section briefly describes the theoretical for- 
malism used to calculate the GT strength distributions 



and the associated electron & positron capture rates. 
The calculated GT± strength distributions are pre- 
sented and compared with measurements and against 
other calculations in Sec. 3. The pn-QRPA calculated 
electron and positron capture rates for iron isotopes 
j|54,55,56pQ^ are presented and explored in Sec. 4. The 
main conclusions of this work are finally presented in 
Sec. 5. 



2 Model Description 

The Hamiltonian for the pn-QRPA calculations was 
taken to be of the form 
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where H'^p is the single-particle Hamiltonian (single 
particle energies and wave functions were calculated in 
the Nilsson model, which takes into account nuclear 
deformations), ^p""" is the pairing force (pairing was 
treated in the BCS approximation), Vqj, is the particle- 
hole (ph) GT force, and Vq^ is the particle-particle 
(pp) GT force. The proton-neutron residual interac- 
tions occurred as particle-hole and particle-particle in- 
teraction. The interactions were given separable form 
and were characterized by two interaction constants x 
and K, respectively. In this work, the values of x ^^nd 
K was taken as 0.15 MeV and 0.07 MeV, respectively. 
Details of choice of these GT strength parameters in 
the pn-QRPA model can be found in Refs. (J32t |33I ). 
Other parameters required for the calculation of weak 
rates are the Nilsson potential parameters, the pairing 
gaps, the deformations, and the Q- values of the reac- 
tions. Nilsson-potential parameters were taken from 



Ref. (|34[ ) and the Nilsson oscillator constant was chosen 
as hijj = 41A~^/^{MeV) (the same for protons and neu- 
trons). The calculated half- lives depend only weakly 
on the values of the pairing gaps (1351 ). Thus, the tradi- 



tional choice of Ap = A„ = 12/y'A{MeV) was applied 
in the present work. The deformation parameter is re- 
cently argued to be one of the most important param- 
eters in pn-QRPA calculations (J36l ) and as such rather 
than using deformations calculated from some theoret- 
ical mass model (as used in earlier calculations of pn- 
QRPA capture rates) the experimentally adopted value 
of the deformation parameters for ^''■^^Fe, extracted by 
relating the measured energy of the first 2+ excited 
state with the quadrupole deformation, was taken from 
Raman et al. (|37t). The incorporation of experimental 
deformations lead to an overall improvement in the cal- 
culation as discussed earlier in Ref. (|28t ). For the case 
of ^^Fe (where measurement lacks) the deformation of 



the nucleus was calculated as 
125(Q2) 
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where Z and A are the atomic and mass numbers, re- 
spectively and Q2 is the electric quadrupole moment 
taken from Ref. (|38l ). Q-values were taken from the 
recent mass compilation of Audi et al. (J39l). 

Capture rates on 54,55,56^g ^q^. ^^^q foUowing two pro- 
cesses mediated by charge weak interaction were calcu- 
lated: 

1. Electron capture 
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2. Positron capture 
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The electron capture (ec) and positron capture (pc) 
rates of a transition from the i*'* state of the parent to 
the j*'* state of the daughter nucleus are given by 
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The value of D was taken to be 6295s (1401 ). B[jS are 
the sum of reduced transition probabilities of the Fermi 
B(F) and GT transitions B(GT). Whereas for 54,56pg 
phonon transitions contribute, in the case of ^^Fe (odd- 
A case) two kinds of transitions are possible. One are 
the phonon transitions, where the odd quasiparticle 
acts as spectator and the other is the transitions of the 
odd quasiparticle itself. In the later case phonon cor- 
relations were introduced to one-quasiparticle states in 
first-order perturbation (|41l). The //^s are the phase 
space integrals. Details of the calculations of phase 
space integrals and reduced transition probabilities can 
be found in Ref. (J19l ). In Eq. (4) {gA/gv)eff is the ef- 
fective ratio of axial and vector coupling constants and 
takes into account the observed quenching of the GT 
strength (|42|). For this project {gA/gv)eff was taken 
(from Ref. (|43l )) as: 
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with {gA/gv)bare = -1.254 (J44l ). Interestiiigly, Vetterli 
and collaborators ([51) and Ronnqvist et al. (J4|) predicted 
the same quenching factor of 0.6 for the RPA calcula- 
tion in the case of ^"^Fc when comparing their measured 
strengths to RPA calculations. 



The total electron (positron) capture rate per unit 
time per nucleus was then calculated using 
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The summation over all initial and final states was car- 
ried out until satisfactory convergence in the rate calcu- 
lations was achieved. Here Pi is the probability of occu- 
pation of parent excited states and follows the normal 
Boltzmann distribution. The pn-QRPA theory allows 
a microscopic state-by-state calculation of both sums 
present in Eq. (6). This feature of the pn-QRPA model 
greatly increases the reliability of the calculated rates 
in stellar matter where there exists a finite probability 
of occupation of excited states. 

In order to further increase the reliability of the cal- 
culated capture rates experimental data were incorpo- 
rated in the calculation wherever possible. In addi- 
tion to the incorporation of the experimentally adopted 
value of the deformation parameter, the calculated exci- 
tation energies (along with their log ft values) were re- 
placed with an experimental one when they were within 
0.5 MeV of each other. Missing measured states were 
inserted and inverse and mirror transitions (if available) 
were also taken into account. No theoretical levels were 
replaced with the experimental ones beyond the exci- 
tation energy for which experimental compilations had 
no definite spin and/or parity. A state-by-state calcula- 
tion of GT± strength was performed for a total of 246 
parent excited states in ^''Fe, 297 states for ^^Fe and 
266 states for ^^Fe. For each parent excited state, tran- 
sitions were calculated for 150 daughter excited states 
using the pn-QRPA model. The band widths of energy 
states were chosen according to the density of states 
to cover an excitation energy of (15-20) MeV in parent 
and daughter nuclei. The summation in Eq. (6) was 
done to ensure satisfactory convergence. The use of a 
separable interaction assisted in the incorporation of a 
luxurious model space of up to 7 major oscillator shells 
which in turn made possible to consider these many 
excited states both in parent and daughter nuclei. 



3 GT-i- strength distributions 

The isovector response of nuclei may be studied us- 
ing the nucleon charge-exchange reactions (p, n) or 
{n,p)] by other reactions such as (•^He,t), (d,^He) or 
through heavy ion reactions. The O" GT cross sec- 
tions (AT = 1,A5 = l,Ai = 0, Ofuv excitations) 
are proportional to the analogous beta-decay strengths. 
Charge-exchange reactions at small momentum transfer 
can therefore be used to study beta-decay strength dis- 
tributions when beta-decay is not energetically possible. 



Table 1 Comparison of measured total GT± strengths 
with microscopic calculations of pn-QRPA and large scale 
shell model in ^'^'''^Fe. For references see text. 
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Experiment 
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7.56 

7.5±0.7 

7.11 



4.26 

3.5±0.7 

3.56 



EGT- 



EGr+ 



10.74 

9.9±2.4 

9.80 



3.71 

2.9±0.3 

2.70 



The {p, n) reactions probes the GT_ strength (corre- 
sponding to beta-minus decay) and the {n,p) reactions 
gives the strength for /3+-decay, i.e. GT_|_ strength. 
The study of (p, n) reactions has the advantage over 
/3-decay measurements in that the GT_ strength can 
be investigated over a large region of excitation energy 
in the residual nucleus. On the other hand the (n,p) 
reactions populates only T = Tq -I- 1 states in all nuclei 
heavier than "^He. This means that other final states 
(including the isobaric analog resonance) are forbidden 
and GT+ transitions can be observed relatively free of 
background. The study of these reactions suggest that 
a reduction in the amount of GT strength is observed 
relative to theoretical calculations. The GT quenching 
is on the order of 30-40 % H). 

In a sense both /3-decay and capture rates are very 
sensitive to the location of the GT_|_ centroid. An 
{n,p) experiment on a nucleus {Z,A) shows where 
in {Z — 1,A) the GT+ centroid corresponding to the 
ground state of {Z,A) resides. Each excited state of 
{Z,A) has its own GT+ centroid in {Z — 1,A) and all 
of these resonances must be included in the stellar rates. 
We do not have the ability to measure these resonances. 
Turning to theory we see that the pioneer calculation 
done by FFN (jlOl ) had to revert to approximations in 
the form of Brink's hypothesis and " back resonances" to 
include all resonances in their calculation. Brink's hy- 
pothesis states that GT strength distribution on excited 
states is identical to that from ground state, shifted 
only by the excitation energy of the state. GT back res- 
onances are the states reached by the strong GT transi- 
tions in the inverse process (electron capture) built on 
ground and excited states. Even the microscopic large- 
scale shell model calculations (jlSJ) had to use the Brink 
assumption to include all states and resonances. On the 
other hand the pn-QRPA model provides a microscopic 
way of calculating the GT_(. centroid and the total GT_|- 
strength for all parent excited states and can lead to a 
fairly reliable estimate of the total stellar rates. 

As a starting point, the model was tested for the case 
of ^''^•^^Fe where measurements of total GT± strength 
are available. As mentioned earlier an overall quench- 
ing factor of 0.6 (|43l ) was adopted in the calculation of 
GT strength in both directions for all isotopes of iron. 
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Fig. 1 Gamow- Teller (GT+ ) strength distributions for 
'"'Fe. The upper panel shows the pn-QRPA results of GT 
strength for the ground state. The middle and lower panels 
show the results for the measured strength distribution (|j) 
and large-scale shell model calculations (jjj), respectively. 
Ej represents daughter states in ^''Mn 



There is a considerable amount of uncertainties present 
in all calculations of stellar rates. The associated un- 
certainties in the pn-QRPA model was discussed in Ref. 
(271). Keeping in mind the uncertainties present also in 
measurements where various energy cutoffs are used as 
a reasonable upper limit on the energy at which GT 
strength could be reliably related to measured AL = 
cross-sections, Table [T] presents the comparison of the 
total GT± strengths against measurements and large- 
scale shell model calculations (J13l ) (referred to as LSSM 
throughout this and proceeding sections) for the case of 
54,56pj,^ Throughout this paper all energies are given in 
units of MeV. For the case of ^^Fe the value of the total 
GT_|_ was taken from Ref. (jj) whereas the total mea- 
sured GT_ was taken from Ref. (jl). The total GT_ 
strength calculated using the pn-QRPA model matches 
very well with the measured strength (Vetterli and col- 
laborators (jsl) also reported a value of 7.5 ± 1.2 as the 
total measured GT_ strength for ^"^Fc). The calculated 
total GT_|_ strength lies close to the upper bound of the 
measured value and are higher than the corresponding 



LSSM calculated strength. For the case of ^^Fe, the 
measured values of the total GT+ was taken from the 
latest measurement by El-Kateb and collaborators (J6|). 
The total measured GT_ was taken from Ref. (J7|). In 
this case the calculated strengths are relatively more 
enhanced as compared to LSSM calculated strengths. 

As mentioned earlier the calculated rates are sen- 
sitive to the location of the GT+ ccntroid. For the 



case of ^''Fe, the calculated ccntroid of 4.06 MeV (|28|) 
is much higher than the LSSM ccntroid of 3.78 MeV. 
The numbers are to be compared to the experimental 
value of 3.7 ± 0.2 MeV which was calculated from the 
measured data presented in Ref. (jj). For ^^Fe, the pn- 
QRPA model calculated the centroid at an excitation 
energy of 3.13 MeV (J28l ) in the daughter nucleus ^^Mn 
whereas LSSM calculated it at an excitation energy of 
2.60 MeV. The centroid extracted from the experimen- 
tal data of El-Kateb et al. (i6|) comes out to be 2.9 ± 
0.2 MeV. Experimental (p,n) data are also available for 
54,56pg^ According to Anderson and collaborators (|8|) a 
large uncertainty exists in the (p,n) measurements and 
the authors were able to perform measurements on ^^Fe 
at 135 MeV with significantly better energy resolution 
than the earlier measurements. Further whereas the 
GT strength in the discrete peaks were extracted accu- 
rately, it was reported that the amount of GT strength 
in the background and continuum remained highly un- 
certain. For the sake of completeness a comparison of 
GT_ centroids is also discussed below. The ccntroid 
of the data of discrete peaks (presented in Table I of 
Ref. fsl)) was calculated to be 7.63 for '^^Fe which is 
to be compared with the much lower pn-QRPA calcu- 
lated value of 5.08. For the case of ^^Fc, the calculated 
centroid from the data reported by Rapaport et al. (|7|) 
comes out to be 8.27. The corresponding pn-QRPA cal- 
culated value is 5.61. Despite greater experimental cer- 
tainties in the (p,n) data one notes that the pn-QRPA 
calculates the centroid at much lower energies in daugh- 
ter nuclei as compared to measurements. However it is 
to be noted that the total strengths (in both direction) 
compares very well with the measured values. Also the 
location of GT-)_ centroid is in very good agreement with 
the measured data which is mainly responsible for the 
calculation of electron capture rates on iron isotopes. 
For reasons not known to the author, LSSM did not 
present the location of GT_ centroids in their paper. 

The GT+ strength distribution for ^"'Fc and ^^Fe are 
depicted in Fig. [1] and Fig. [21 respectively. The up- 
per panel displays the calculated GT+ distribution us- 
ing the pn-QRPA model. The middle panel shows the 
measured data and the bottom panel gives the corre- 
sponding LSSM strength distributions. One notes that 
an appreciable measured strength lies in the daughter 
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Fig. 2 Gamow- Teller (GT+ ) strength distributions for 
^^Fe. The upper panel shows the pn-QRPA results of GT 
strength for the ground state. The middle and lower panels 
show the results for the measured values (|a) and large-scale 
shell model calculations (jjj), respectively. Ej represents 
daughter states in ^®Mn. 

excitation energies up to around 10 MeV. The LSSM 
strength peaks at much lower excitation energy. The 
pn-QRPA distribution spreads over higher excitation 
energies akin to measured strength with a peak around 
7.1 MeV for the case of ^^Fe. For the case of ^^^Fe 
the pn-QRPA calculated spread is not as good as in 
case of ^^Fe. Nonetheless the pn-QRPA calculates its 
peak around 4 MeV much higher in energy than the 
LSSM peak accumulating all the measured higher lying 
strength in a narrow resonance region. 

Fig. [3] shows the calculated GT+ strength distribu- 
tion for ^^Fc. Here no experimental GT distribution 
was available for comparison. Instead the calculated 
strength for ground and first two excited states of '"'^Fe 
(at 0.41 MeV and 0.93 MeV, respectively) is displayed. 
It may be seen from Fig. 3 that the GT strength is frag- 
mented over many daughter excited states and peaks 
at relatively high excitation energies (around 8 MeV) 
in the daughter ^^Mn. 

Figs, m [S] and [5] display the GT_ strength distri- 
butions for the ground and first two calculated excited 



Fig. 3 Gamow- Teller (GT+ ) strength distributions for 
^^Fe. The upper panel shows the pn-QRPA results of GT 
strength for the ground state. The middle and lower panels 
show the results for the calculated first and second excited 
states of ^^Fe. Ej represents daughter states in ^^Mn. 



states of ^^Fe, ^^Fe and ^^Fe, respectively. One notes 
the enhancement in the total strength in the GT_ di- 
rection. For the case of ^^Fe the peaks shift to higher 
excitation energy in daughter for parent excited states. 
Correspondingly the energy centroids shift to higher en- 
ergy for these excited states. For the case of ^^Fe one 
notes that bulk of the strength resides in a narrow reso- 
nance region around 10 MeV in daughter. In the case of 
^^Fc the low-lying peaks become considerably small in 
magnitude for the parent excited states. Figs.[Tl[2l[3l|4l 
El and [6] confirm that the calculated strength is frag- 
mented and well spread. At low temperatures and den- 
sities these low-lying discrete strengths may very well 
dominate the rates and play an important role in the 
state-by-state evaluation of both sums (see Eq. (6)). 

The ground and excited states strength functions cal- 
culated within the framework of the pn-QRPA model 
are required in the microscopic calculation of stellar 
capture rates. The ASCII files of the GT± strength 
distributions for all parent excited states are available 
and can be requested from the author. 
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Fig. 4 Gamow- Teller (GT-) strength distributions for 
^''Fe. The upper panel shows the pn-QRPA results of GT 
strength for the ground state. The middle and lower panels 
show the results for the calculated first and second excited 
states of ^^Fe. Ej represents daughter states in ^''Co. 



Fig. 5 Gamow- Teller (GT_) strength distributions for 
^^Fe. The upper panel shows the pn-QRPA results of GT 
strength for the ground state. The middle and lower panels 
show the results for the calculated first and second excited 
states of ^^Fe. Ej represents daughter states in ^^Co. 



4 Electron and positron capture rates 

This section presents and explores the results of the cal- 
culated electron and positron capture rates on 54,55,56pg 
in stellar environment and also compare the pn-QRPA 
capture rates with the LSSM calculation and against 
the pioneering calculation of FFN (llOl ). As discussed 
earlier ground and excited states GT strengths calcu- 
lated earlier contribute in the state-by-statc calculation 
of electron and positron capture rates (via Eq. 4) . 

Figs. [3 |8] and |9] show the calculated electron cap- 
ture rates on ^^'^^•^^'Pe. These figures show four panels 
depicting the calculated electron capture rates at se- 
lected temperature and density domain. The upper left 
panel shows the electron capture rates in low-density re- 
gion (p[gcm^^] = lO"-^, 10^-^ and 10^-^), the upper right 
in medium-low density region {p[gcm^'^] = lO'^'^, 10*'^ 
and 10^'^), the lower left in medium-high density re- 
gion {p[gcm~^] = 10^-^, 10'^-^ and 10*-^) and finally the 
lower right panel depicts the calculated electron capture 
rates in high density region {p[gcm~^] ~ 10^'^, 10^°'^ 
and 10^^). The electron capture rates are given in loga- 



rithmic scales (to base 10) in units of s^^. Tg gives the 
stellar temperature in units of 10^ K. Figs.[7]and IHlare 
similar in nature depicting the electron capture rates on 
even-even isotopes of iron. The rates are relatively en- 
hanced for the case of ^^Fe by orders of magnitude (for 
the first three panels). It can be seen from these figures 
that in the low density regions the capture rates, as a 
function of temperature, arc more or less superimposed 
on one another. This means that there is no appre- 
ciable change in the rates when increasing the density 
by an order of magnitude. However as one goes from 
the medium-low density region to high density region 
these rates start to 'peel off' from one another. Or- 
ders of magnitude difference in rates are observed (as a 
function of density) in high density region. For a given 
density the rates increase monotonically with increas- 
ing temperatures. For all three isotopes of iron, the 
calculated electron capture rates are noticeably bigger 
than the competing /S"*" decay rates and dominate the 
weak rates for charge-decreasing nuclear transitions in 
stellar matter. 
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Fig. 6 Gamow- Teller (GT-) strength distributions for 
®^Fe. The upper panel shows the pn-QRPA results of GT 
strength for the ground state. The middle and lower panels 
show the results for the calculated first and second excited 
states of ^^Fe. Ej represents daughter states in ^®Co. 



The effects of positron capture are estimated to be 
smaller for stars in the mass range (10 < Mp, < 40) 



and could be greater for more massive stars (|3lf). The 
positron capture rates arc normally bigger than the 
competing /3~ rates for all isotopes of iron specially 
at high temperatures. Figs, [ini HH and [T^] again 
show four panels depicting the calculated positron cap- 
ture rates at selected temperature and density do- 
main. The upper left panel shows the positron cap- 
ture rates in low-density region {p[gcm~^] = 10°'^, 10^'^ 
and 10^'^) ; the upper right in medium-low density re- 
gion {p[gcm-^] = 103■^10''■^ and lO^-^), the lower left 
in medium- high density region {p[gcm^^] = 10^'^, 10^'^ 
and 10^'^) and finally the lower right panel depicts the 
calculated positron capture rates in high density region 
{p[gcm-^] = 109■^10l°■^ and 10"). The positron cap- 
ture rates are given in logarithmic scales in units of 
s~^. Tg gives the stellar temperature in units of 10^ 
K. One should note the order of magnitude enhance- 
ment in positron capture rates as the stellar tempera- 
ture increases. Around presupernova temperatures the 
positron capture rates are very slow as compared to 
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Fig. 7 (Color online) Electron capture rates on Fe, as a 
function of stellar temperatures, for different selected den- 
sities. Temperatures are given in 10^ K. Densities are given 
in units of gcmT and logioAec represents the log of electron 
capture rates in units of sec^^ . 

the corresponding electron capture rates. When the 
temperature of the stellar core increases further the 
positron capture rates shoot up. The positron capture 
rates decrease with increasing densities, in contrast to 
the electron capture rates which increase as density in- 
creases. As temperature rises or density lowers (the de- 
generacy parameter is negative for positrons) , more and 
more high-energy positrons are created leading in turn 
to higher capture rates. It is worth mentioning that 
the positron capture rates are very small numbers and 
can change by orders of magnitude by a mere change 
of 0.5 MeV, or less, in parent or daughter excitation 
energies and are more reflective of the uncertainties in 
the calculation of energies. 

In order to present an estimate of the contribution of 
excited states to the total capture rate, the ground state 
capture rate and the ratio of this rate to the total rate 
were computed at selected points of stellar temperature 
and density. These contributions vary for different iso- 
topes and lead to interesting consequences as presented 
below. 

Tables El El and |4] show the contribution of ex- 
cited states in the calculation of total capture rates on 
54,55,56pg^ In all tables the capture rates and the ratios 
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are calculated at selected points of density and tem- 
perature shown in first column. Within the parenthesis 
the first number gives the density in units of gcmT^ and 
the second number denotes the stellar temperature in 
units of 10^ K. In the second column Aec(G) gives the 
ground state electron capture rates whereas Rec{G/T) 
denotes the ratios of the ground-state electron capture 
rate to the total rate. The ground state positron cap- 
ture rates, Apc(G), and the ratios of the ground-state 
positron capture rate to the total rate, Rpc{G/T), are 
given in the fourth and fifth column, respectively. All 
capture rates are given in units of sec~^. 

The results for ^"^Fe are very interesting and are de- 
picted in Table [H The total electron capture rate in- 
creases with increasing temperatures and densities. At 
low temperatures and densities (Tg = 1 and pYe < 
lO^gcTO"^) one notes that almost all of the contribu- 
tion to the total electron capture rate comes from the 
excited states. At a first glance this might look odd. 
To explain this result one has to go back to Eq. (4) 
which tells that the partial capture rate is a product of 
three factors: a constant, phase space integrals (which 
are functions of temperature and density) and reduced 
transition probabilities. For a particular parent state 
(i) all such transitions are summed over daughter states 
(j) and then multiplied by the probability of occupa- 



tion of that parent excited state. These partial capture 
rates are finally summed over all parent excited states 
to calculate the total capture rate as a function of stel- 
lar temperature and density (Eq. (6)). For the given 
physical conditions, the ground state electron capture 
rate ^ 10~^*s~^. Probability of occupation of ground 
state is essentially 1 at such low temperatures. However 
for the first excited state (at 1.41 MeV) the partial rate 
increases roughly by 10 orders of magnitude to around 
10~®s~^. This increase is attributed to a simultane- 
ous increase in phase space (= Q + Ei — Ej where Q 
= -0.697 MeV (|39|)) as well as an increase in the to- 
tal GT+ strength from 4.26 units to 5.12 units (refer 
to Table 2 of Ref. ([H)). Even after multiplying by 
the much smaller probability of occupation of parent 
excited state (^ lO^'*) the partial capture rate from 
the first excited state is around 3 orders of magnitude 
bigger than that from ground state at low tempera- 
tures. As density increases the phase space from excited 
states decreases resulting in a larger contribution to to- 
tal rate from ground state. In the limiting case of densi- 
ties around 10^^ gcm~^ all the contribution to the total 
rate comes from ground state. At higher temperatures 
(3 < Tg < 10) and densities (< lO'^gcm^^) the contri- 
bution of ground state to total rate is roughly 50%. At 
still higher temperatures (Tg = 30) the ground state 
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Fig. 10 (Color online) Same as Fig. [7] but for positron Fig. 11 (Color online) Same as Fig. [7| but for positron 
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contributes around 6-7% at all densities. At Tg = 30, 
the partial capture rates are significant also for higher 
lying excited states (as probability of occupation of 
these high-lying states then becomes appreciably high) 
and contribute effectively to the total capture rate. For 
the corresponding contribution of excited states to to- 
tal positron capture rate, one sees that these contribu- 
tions arc essentially dependent on stellar temperatures 
and independent of stellar densities. Further one notes 
that the contribution of excited states dominate in total 
positron capture rates. These features are reflective of 
the fact that the positron capture rates are dominated 
by phase space integrals. For a given temperature the 
total positron capture rate increases by a constant fac- 
tor of the corresponding ground state rate independent 
of the densities. At low temperatures the ground state 
contribution is almost negligible for reasons mentioned 
above. As the stellar temperature increases so does the 
ground state contribution until it reaches a maximum 
around Tg = 10 of the order of 27%. At Tg = 30 
roughly 95% of the contribution comes from excited 
states. 

For the case of ^^Fe the story is much different (Ta- 
ble [3]). Even though the first excited state is at 0.41 
MeV (as against 1.41 McV for the case of ^"'Fc) resulting 
in a much bigger probability of occupation, the partial 



electron capture rate from first excited state is around 
7 orders of magnitude smaller than the corresponding 
ground state electron capture rate ~ 10~^''s~^ (the to- 
tal GT+ strength decreases for the first excited state 
from 4.68 units to 4.43 units (|28l )). As a result the total 
electron capture rate is determined almost entirely from 
ground state rate at low temperatures. This is in sharp 
contrast to the previous case. The low-lying excited 
states do play their role at higher temperatures and 
densities. The behavior of ground state contribution 
to the total electron capture rate is as per expectation 
and decreases appreciably as the stellar temperature 
increases. Again at Tg = 30 almost 97% of the contri- 
bution is from the parent excited states at all densities 
akin to the case of ^''Fe. The excited states contribu- 
tion to the total positron capture rate follows a similar 
trend as in previous case with a maximum contribu- 
tion of around 25% from ground state at Tg = 10. At 
Tg = 3 the ground state contributes ~ 2% (as against 
17% in the case of ^^Fe). At much higher temperature 
(Tg ~ 30) this contribution is again around 2%. 

The excited states contribution to total electron 
capture rate on ^^Fe at low temperatures and densi- 
ties falls roughly in between the two extreme cases of 
^''Fc and ^SFe (Table S]). Here roughly 50% contri- 
bution comes from the excited states at low temper- 



12 








y"^^^^ 


,^ 


-15 




II) 
















r^ 


-30 


T 








0) 

n 




1 — p=10°-= 




-4b 


1 -.-p=io^-= - 




-60 


- ^ p=10" 



5 10 15 20 25 30 






^ ^ " 




/"^ 


^-^ -15 


/ 


(D 


f 


a. 



45 




-.-P=10'= 
-.-p=10^-=" 


-60 




^P=10== 





-20 

91 

-i -40- 








5 10 15 20 25 30 




W) 



Fig. 12 (Color online) Same as Fig. [7] but for positron 
capture rates on ^^Fe. 



atures and densities. Example giving at Tg = 1 and 
pYg < lQgcm~^, the total ground state electron cap- 
ture rate is around 10~^^s~^. The contribution from 
first excited state (0.85 MeV) is around lO^^^g-i ^his 
is again attributed partly to the increased total GT+ 
strength of 5.15 units as against 3.71 units (for ground 
state) (|28l ). This rate when multiplied with the occu- 
pation probability of around 10~'^ gives a partial elec- 
tron capture rate of roughly equal magnitude as that 
from ground state. The excited states contribution in- 
creases at higher temperatures with roughly 95% of 
the contribution coming from parent excited states at 
Tg = 30. One also notes that the ground state contri- 
bution to the total electron capture rates increases at 
pYe ~ lO^^gcm^^. Specially at low temperatures it is 
almost totally dictated by ground state rates. This is 
because ground and excited state rates are of compa- 
rable magnitudes at this density. However the proba- 
bility of occupation of excited states is smaller by or- 
ders of magnitude at lower temperatures. The excited 
states contribution to total positron capture rate fol- 
lows a similar trend as in previous cases with a maxi- 
mum contribution of around 30% from ground state at 
Tg = 10. Tables H El and H highlight the contribution 
of partial capture rates from parent excited states to the 
total capture rate and show disparate behavior of these 



Table 2 The ground state electron and positron capture 
rates, Xec{G), Xpc{G) for ^*_Fe in units of sec^^ . Given 
also are the ratios of the ground state capture rates to total 
capture rates, Rec{G/T), Rpc(G/T). The first column gives 
the corresponding values of stellar density, pYe {gcm~^), 
and temperature, Tg (in units of 10^ K), respectively. 







54 Fe 






(pYe.Tg) 


Aec(G) 


Rec(G/T) 


Apc(G) 


Rpc(G/T) 


(10,1) 


2.89E-18 


1.61E-03 


3.32E-51 


2.83E-04 


(10,3) 


7.49E-08 


4.10E-01 


5.44E-18 


1.71E-01 


(10,10) 


1.16E-02 


4.69E-01 


1.53E-05 


2.69E-01 


(10,30) 


2.89E+00 


5.77E-02 


3.03E-01 


5.22E-02 


(10^,1) 


4.88E-18 


1.61E-03 


1.97E-51 


2.83E-04 


(10^,3) 


7.51E-08 


4.10E-01 


5.44E-18 


1.71E-01 


(10^,10) 


1.16E-02 


4.69E-01 


1.53E-05 


2.69E-01 


(10^,30) 


2.89E+00 


5.77E-02 


3.03E-01 


5.22E-02 


(10', 1) 


3.22E-12 


7.91E-02 


2.97E-57 


2.83E-04 


(10'',3) 


3.90E-06 


4.96E-01 


1.05E-19 


1.71E-01 


(10^,10) 


1.45E-02 


4.71E-01 


1.22E-05 


2.69E-01 


(10', 30) 


2.92E+00 


5.77E-02 


3.01E-01 


5.22E-02 


(10^\l) 


2.56E+04 


l.OOE+00 


7.98E-172 


2.83E-04 


(10^1,3) 


2.56E+04 


9.95E-01 


3.49E-58 


1.71E-01 


(10",10) 


1.93E-f04 


7.12E-01 


1.49E-17 


2.69E-01 


(10",30) 


3.14E+03 


7.06E-02 


4.12E-05 


5.20E-02 



Table 3 f 


)ame as Table [2] but for 


"^Fe. 




(pYe,T9) 


Aec(G) 


""be 

Rec(G/T) 


Apo(G) 


Rpc(G/T) 


(10,1) 


3.85E-10 


9.98E-01 


1.71E-32 


2.48E-08 


(10,3) 


6.85E-07 


3.56E-01 


1.47E-12 


2.01E-02 


(10,10) 


2.93E-03 


1.53E-01 


1.72E-04 


2.35E-01 


(10,30) 


1.51E+00 


2.56E-02 


3.61E-01 


2.34E-02 


(10'',1) 


6.51E-10 


9.98E-01 


l.OlE-32 


2.48E-08 


(10^,3) 


6.86E-07 


3.56E-01 


1.47E-12 


2.01E-02 


(10''', 10) 


2.93E-03 


1.53E-01 


1.72E-04 


2.35E-01 


(10-', 30) 


1.51E+00 


2.56E-02 


3.62E-01 


2.34E-02 


(10',1) 


1.22E-05 


9.99E-01 


1.53E-38 


2.48E-08 


(10'',3) 


2.90E-05 


4.28E-01 


2.83E-14 


2.01E-02 


(10'',10) 


3.68E-03 


1.54E-01 


1.37E-04 


2.35E-01 


(10''',30) 


1.52E+00 


2.56E-02 


3.59E-01 


2.34E-02 


(10^\l) 


2.25E+04 


9.92E-01 


4.11E-153 


2.48E-08 


(10^^,3) 


1.83E+04 


8.07E-01 


9.44E-53 


2.01E-02 


(10",10) 


8.43E+03 


3.35E-01 


1.67E-16 


2.35E-01 


(10",30) 


2.08E+03 


4.03E-02 


4.94E-05 


2.33E-02 



contributions for all three cases. This analysis further 
stresses on the fact that a microscopic calculation of GT 
strength function from excited states is indeed required 
for a reliable estimate of the total capture rates. 

One relevant question would then be to ask how the 
calculated rates compare with other prominent calcula- 
tions. For the comparison two important calculations of 
capture rates were taken into consideration: the pioneer 
calculation of FFN (|lO( ) which is still used in many sim- 
ulation codes and the microscopic calculation of LSSM. 

The comparisons are presented in a tabular form. 
Tables [51 [HI and [71 show the comparison of calculated 
electron capture rates with those of FFN and LSSM for 
54pg^ 55pg g^jj^ ^^Fc, respectively. Here the ratios of the 
calculated electron capture rates to those of FFN and 
LSSM are presented at selected temperature and den- 
sity points. For the case of ^*Fe, the calculated electron 
capture rates are in good agreement with the LSSM re- 
sults specially at high temperatures and densities. 

At low temperatures and densities the reported elec- 
tron capture rates are bigger than the LSSM capture 
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Table 5 Ratios of calculations of electron capture rates on ''"^Fe at different selected densities and temperatures. QRPA 
implies the reported rates whereas SM and FFN denote rates calculated by Ref. (|l3h and Ref. (|lOr). respectively. 
Tg QRPA/SM QRPA/FFN QRPA/SM QRPA/FFN QRPA/SM QRPA/FFN QRPA/SM QRPA/FFN 



lOgcm 



lOgcm 



10' gem 



10 gem 



10 gem 



10 gc77i 



lO^^gcm" 



10 gem" 



1 


4.49E+00 


1.04E-01 


4.49E+00 


1.05E-01 


3.05E+00 


1.31E-01 


8.59E-01 


3.00E-01 


3 


2.50E+00 


1.93E-01 


2.50E+00 


1.93E-01 


2.38E+00 


2.36E-01 


8.59E-01 


2.99E-01 


10 


9.95E-01 


4.36E-01 


9.93E-01 


4.36E-01 


9.95E-01 


4.37E-01 


8.45E-01 


2.99E-01 


30 


1.14E+00 


4.00E-01 


1.14E+00 


4.00E-01 


1.14E+00 


4.00E-01 


1.05E+00 


3.62E-01 



Table 4 S 


)ame as Table [2] but for 


5«Fe. 




(pYe.Ts) 


Aoc(G) 


Rec(G/T) 


Apc(G) 


R.pc(G/T) 


(10,1) 


8.33E-27 


4.70E-01 


6.49E-37 


7.01E-08 


(10,3) 


1.59E-10 


5.69E-01 


2.26E-13 


1.77E-02 


(10,10) 


1.60E-03 


4.07E-01 


3.80E-04 


3.21E-01 


(10,30) 


1.52E+00 


5.33E-02 


9.90E-01 


5.77E-02 


(10^1) 


1.41E-26 


4.70E-01 


3.84E-37 


7.01E-08 


(10^,3) 


1.59E-10 


5.69E-01 


2.25E-13 


1.77E-02 


(10^,10) 


1.60E-03 


4.07E-01 


3.80E-04 


3.21E-01 


(10^,30) 


1.52E+00 


5.33E-02 


9.91E-01 


5.77E-02 


(10', 1) 


9.30E-21 


4.77E-01 


5.82E-43 


7.01E-08 


(10^3) 


8.28E-09 


5.71E-01 


4.34E-15 


1.77E-02 


(io'',io) 


2.01E-03 


4.07E-01 


3.03E-04 


3.20E-01 


(10'',30) 


1.53E+00 


5.33E-02 


9.83E-01 


5.77E-02 


(lo^M) 


1.80E+04 


l.OOE+00 


1.56E-157 


7.01E-08 


(10",3) 


1.74E+04 


9.64E-01 


1.45E-53 


1.77E-02 


(10",10) 


1.14E+04 


6.03E-01 


3.70E-16 


3.18E-01 


(10",30) 


2.13E+03 


6.35E-02 


1.35E-04 


5.74E-02 



rates by around a factor of four (Table [5]) . During the 
oxygen shell burning and silicon core burnin g of mas- 
sive stars (as per the simulation results of (|31l)) the 
pn-QRPA calculated rates are enhanced roughly by a 
factor of three as compared to LSSM rates. Whereas 
the individual discrete transitions between initial and 
final states matter at low temperatures and densities, 
it is the total GT strength that counts at high temper- 
atures and densities. It is again reminded that Brink's 
hypothesis is not assumed in this calculation (which was 
adopted in LSSM calculation of weak rates). (Also see 
the discussion on contribution of excited state partial 
rates above.) The ground state pn-QRPA centroid is 
placed at much higher energy in ^^Mn as compared to 
LSSM centroid (|28l). However, pn-QRPA calculates a 
much bigger total strength. FFN rates are up to around 
an order of magnitude enhanced compared to pn-QRPA 
rates at low temperatures and densities. FFN did not 
take into effect the process of particle emission from 
excited states and their parent excitation energies ex- 
tended well beyond the particle decay channel. These 
high lying excited states begin to show their effect spe- 
cially at high temperatures and densities. Also one 
notes that FFN neglected the quenching of the total 
GT strength. The LSSM electron capture rates were 
smaller than the FFN rates by, on average, an order of 
magnitude (J3ll ). 

The electron capture rates on ^^Fe arc most effective 
during the oxygen shell burning till around the igni- 
tion of the first stage of convective silicon shell burning 
of massive stars (J3lh . For the corresponding temper- 



atures and densities, the pn-QRPA rates are in very 
good agreement with the LSSM rates at all tempera- 
tures and densities (Table |6]) . This excellent agreement 
is attributed to the fact that the total rates for ^^Fe are 
commanded by ground state contribution (see Table [3|) 
and both pn-QRPA and LSSM perform a microscopic 
calculation of the ground-state GT strength function. 
The comparison with FFN is good at low temperatures 
only (Tg ^1). At higher temperatures FFN rates are 
bigger than the LSSM and pn-QRPA rates by an order 
of magnitude for reasons mentioned before. 

The average comparison of calculated electron cap- 
ture rates is again very good against the LSSM rates 
for the case of ^^Fe (Table [7]). However at low temper- 
atures and densities the LSSM electron capture rates is 
around 10% enhanced as compared to pn-QRPA num- 
bers. As mentioned above the comparison of the re- 
ported and LSSM electron capture rates may be traced 
back to the calculation of excited state partial capture 
rates (Tables [2l |3l and |4]). One notes that at Tg = 1 
and pYe < lO^cm"'^ the electron capture rates were 
dominated by excited state contributions for the case of 
^"'Fe. For the case of ^^Fe it was entirely dominated by 
ground state contribution and the case of ^^Fe had a SO- 
SO contribution. Correspondingly the comparison with 
LSSM rates is excellent for the case of ^^Fe. The pn- 
QRPA rates arc enhanced by a factor of 4 for ^^Fe and 
some fluctuations in rate comparison can be seen for 
the case of ^^Fe. The electron capture rates on ^^Fc are 
very important for the pre-supernova phase of massive 
stars. FFN rates are again enhanced by up to an or- 
der of magnitude as compared to pn-QRPA and LSSM 
rates. 

Tables mini and [TUlshow the corresponding compari- 
son for the calculated positron capture rates. Again the 
pn-QRPA calculated positron capture rates are com- 
pared with the LSSM and FFN rates. It is reminded 
that the positron capture rates are smaller than the 
corresponding electron capture rates by orders of mag- 
nitude and these small numbers can change appreciably 
by a mere change of O.S MeV in phase space and are ac- 
tually more reflective of the uncertainties in calculation 
of the energy eigenvalues (for both parent and daughter 
states). Further it is also evident from Tables HI IHl and 
1101 that positron capture rates have a dominant con- 
tribution from parent excited states (the ground state 
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Table 6 Same as Table 5, but for electron capture on Fe. 
To 



QRPA/SM 
lOgcni^' 



QRPA/FFN 
lOgcni^' 



QRPA/SM 



QRPA/FFN QRPA/SM 



QRPA/FFN 



QRPA/SM 



QRPA/FFN 



1 


l.OlE+00 


l.OOE+00 


l.OlE+00 


9.98E-01 


1.02E+00 


l.OlE+00 


9.04E-01 


2.04E-01 


3 


1.45E+00 


3.72E-01 


1.44E+00 


3.72E-01 


1.44E+00 


3.58E-01 


8.91E-01 


2.04E-01 


10 


1.12E+00 


1.12E-01 


1.12E+00 


1.12E-01 


1.12E+00 


1.12E-01 


9.23E-01 


2.19E-01 


30 


1.72E+00 


3.50E-01 


1.72E+00 


3.51E-01 


1.72E+00 


3.50E-01 


1.44E+00 


3.53E-01 



Table 7 Same as Table 5, but for electron capture on Fe. 
Tg QRPA/SM 



lOgcm 



QRPA/FFN 
lOgcm^' 



QRPA/SM 
W^gcm^^ 



QRPA/FFN QRPA/SM 
lO^gcm"^ lO'^gcm"^ 



QRPA/FFN 
lO'^gcm"-'^ 



QRPA/SM QRPA/FFN 
lO"^"^ gem" 



-,11 



10 gem'" 



1 

3 

10 
30 



8.87E-01 
1.12E+00 
1.02E+00 
1.29E+00 



2.77E-01 
2.63E-01 
4.80E-01 
5.04E-01 



8.87E-01 
1.12E+00 
1.02E+00 
1.29E+00 



2.78E-01 
2.63E-01 
4.80E-01 
5.04E-01 



l.OlE+00 
1.17E+00 
1.02E+00 
1.29E+00 



3.37E-01 
2.77E-01 
4.80E-01 
5.02E-01 



1.13E+00 
l.lOE+00 
1.07E+00 
1.37E+00 



3.37E-01 
3.33E-01 
3.25E-01 
4.14E-01 



contributes at the maximuni by one third). Looking 
at Tabic [5] for the case of ^^Fe one notes that the pn- 
QRPA calculated positron capture rates are suppressed 
by up to 5 orders of magnitude as compared to LSSM 
and FFN rates at low temperatures and densities. The 
comparison improves as the temperature and density 
increases. At Tg = 30 the reported rates are in fact 
enhanced by a factor of 8 as compared to LSSM calcu- 
lated rates and are in reasonable agreement with FFN 
rates. 

The comparison follows a similar trend for the case 
of 55pg (Tabled. Here the pn-QRPA calculated rates 
are suppressed by as much as 3 orders of magnitude 
compared with the other calculations at low tempera- 
tures. The comparison improves as the density of stellar 
core stiffens. The reported rates arc enhanced at higher 
temperatures by as much as an order of magnitude as 
compared to LSSM numbers. FFN rates are again in 
reasonable agreement with reported rates at Tg = 30. 

The comparison with LSSM rates improves as one 
matches the results from ^^Fe to ^^Fe and finally to 
^^Fe. In the case of ^^Fe the reported positron capture 
rates are suppressed by around 2 orders of magnitude 
at Tg = 1 (Table [TU)) . At higher densities and tempera- 
tures the rates are in very good comparison. The FFN 
rates are enhanced by roughly 4 orders of magnitude at 
low temperatures and densities. The comparison im- 
proves as the density and temperature of stellar core 
increases. 



5 Summary and conclusions 

The microscopic calculation of Gamow- Teller strength 
distributions (GT±) of three key isotopes of iron of as- 
trophysical importance was presented. The calculated 
strengths were in good comparison to the measured 
strengths (for the case of ^^'^^Fe). The pn-QRPA calcu- 
lated total GT strengths were greater in magnitude as 



compared to those using LSSM. The results also high- 
lighted the fact that the Brink's hypothesis and back 
resonances may not be a good approximation to use in 
stellar calculation of weak rates. 

The pn-QRPA model was also used to calculate the 
associated electron and positron capture rates of these 
isotopes of iron with astrophysical importance. Defor- 
mations of nuclei were taken into account and for the 
first time the reported calculation took into considera- 
tion the experimental deformations. The rates are cal- 
culated on an extensive temperature-density grid point 
suitable for interpolation purposes that might be re- 
quired in collapse simulations. The electronic versions 
of these files may be requested from the author. 

During the oxygen and silicon core burning phase 
of massive stars the pn-QRPA electron capture rates 
on ^"'Fc are around three times bigger than those calcu- 
lated by LSSM. The comparison with LSSM gets better 
for proceeding pre-supernova and supernova phases of 
stars. The pn-QRPA calculated electron capture rates 
on 55,56pg g^j.g jj-^ overall excellent agreement with the 
LSSM rates. The calculated positron capture rates are 
generally smaller (by as much as five orders of magni- 
tude) at low temperatures and densities. The calcula- 
tion further discourages any noticeable contribution of 
positron capture rates on iron isotopes for the physics 
of core-collapse. 

Due to the weak interaction processes (capture and 
decay rates) the value of Icpton-to-baryon (Ye) ratio for 
a massive star changes from 1 (during hydrogen burn- 
ing) to roughly 0.5 (at the beginning of carbon burn- 
ing) and finally to around 0.42 just before the collapse 
to a supernova explosion. The temporal variation of Yg 
within the core of a massive star has a pivotal role to 
play in the stellar evolution and a fine-tuning of this pa- 
rameter at various stages of presupernova evolution is 
the key to generate an explosion. The electron capture 
tends to reduce this ratio whereas the positron capture 
increases it. This paper reported on the microscopic 
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Table 8 Same as Table 5, but for positron capture on ^ Fe. 



Ta 



QRPA/SM 
lOgcni^' 



QRPA/FFN QRPA/SM 



QRPA/FFN QRPA/SM 



QRPA/FFN QRPA/SM QRPA/FFN 



10 gem 



10 gem" 



10 gcm~ 



1 


3.50E-05 


3.53E-05 


3.50E-05 


3.53E-05 


3.50E-05 


3.52E-05 


9.98E-01 


9.98E-01 


3 


2.91E-02 


3.18E-02 


2.91E-02 


3.17E-02 


2.91E-02 


3.18E-02 


2.90E-02 


3.18E-02 


10 


6.97E-01 


8.47E-01 


6.97E-01 


8.47E-01 


6.97E-01 


8.43E-01 


6.89E-01 


8.24E-01 


30 


7.59E+00 


8.51E-Q1 


7.60E+00 


8.53E-01 


7.62E+00 


8.55E-01 


7.41E+00 


8.38E-01 



Table 9 Same as Table 5, but for positron capture on '^^ Fe. 
Tg QRPA/SM 



lOgcm 



QRPA/FFN 
lOgcm^' 



QRPA/SM 
10'''gcm"''' 



QRPA/FFN QRPA/SM 
lO^gcm"^ lO^'gcm"'' 



QRPA/FFN QRPA/SM QRPA/FFN 



10 gcva 



10 gem" 



10 gem" 



1 


1.85E-03 


1.57E-03 


1.85E-03 


1.57E-03 


1.85E-03 


1.57E-03 


9.98E-01 


9.98E-01 


3 


4.73E-01 


4.12E-02 


4.73E-01 


4.12E-02 


4.72E-01 


4.06E-02 


4.69E-01 


4.06E-02 


10 


1.45E+00 


1.34E-01 


1.45E+00 


1.34E-01 


1.44E+00 


1.34E-01 


1.40E+00 


1.31E-01 


30 


1.15E+01 


8.39E-01 


1.15E+01 


8.41E-01 


1.15E+01 


8.39E-01 


1.12E+01 


8.24E-01 



calculation of electron and positron capture rates on 
iron isotopes and also highlighted the differences with 
previous key calculations. What affect the calculated 
rates may have on the simulation result is difficult to 
predict at this stage. Core-collapse simulators are urged 
to check the affect of incorporating pn-QRPA rates in 
their simulation codes for possible interesting outcome. 
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Table 10 Same as Table 5, but for positron capture on ^®Fe. 



Tg QRPA/SM 



QRPA/FFN 
lOgcni^' 



QRPA/SM 



QRPA/FFN 
10^ gem"'' 



QRPA/SM 
10' gcm^' 



QRPA/FFN 



QRPA/SM 
lO^^gcm^^ 



QRPA/FFN 
10"ffcm"^ 



1 


4.12E-02 


5.32E-04 


4.12E-02 


5.32E-04 


4.11E-02 


5.32E-04 


9.98E-01 


9.98E-01 


3 


2.36E+00 


7.05E-02 


2.36E+00 


7.06E-02 


2.35E+00 


7.03E-02 


2.33E+00 


7.03E-02 


10 


1.84E+00 


1.07E-01 


1.84E+00 


1.07E-01 


1.84E+00 


1.06E-01 


1.81E+00 


1.03E-01 


30 


8.15E+00 


4.79E-01 


8.17E+00 


4.80E-01 


8.15E+00 


4.79E-01 


7.96E+00 


4.67E-01 
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